Directed transport born from chaos in asymmetric antidot structures 
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It is shown that a polarized microwave radiation creates directed transport in an asymmetric 
antidot superlattice in a two dimensional electron gas. A numerical method is developed that allows 
to establish the dependence of this ratchet effect on several parameters relevant for real experimental 
studies. It is applied to the concrete case of a semidisk Galton board where the electron dynamics 
is chaotic in the absence of microwave driving. The obtained results show that high currents can 
be reached at a relatively low microwave power. This effect opens new possibilities for microwave 
control of transport in asymmetric superlattices. 
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INTRODUCTION 

The appearance of a directed transport induced by ra- 
diation in asymmetric systems is known as the photogal- 
vanic effect. By this effect the radiation creates charge 
transport in the bulk of the asymmetric structure in ab- 
sence of any applied cic-voltage. The theoretical investi- 
gations of this phenomena have been started almost 30 
years ago in Refs. 0, 0]- The interest to this subject 
has been renewed recently with the studies of ratchets 
that appear when a system is displaced from thermal 
equilibrium by a periodic variation of system parameters 
(for reviews see Refs. 0, 0). One of surprising prop- 
erties of ratchets is that directed transport can emerge 
in presence of a zero mean force. This phenomenon has 
a generic origin and appears in various physical systems 
including vortices in Josephson junction arrays [EJ(J0> 
cold atoms |8|, ma crop orous silicon membranes |9(, mi- 
crofluidic channels [lfj and other systems. 

Nowadays technology allows to prepare artificial anti- 
dot superlattices in semiconductor heterostructure with 
two dimensional electron gas (2DES). The conduction 
properties of these samples has been tested in experi- 
ments [TlL ll'J which showed an important contribution 
of periodic orbits. The structure of these superlatticies 
is similar to a periodic lattice of rigid disks placed on a 
plane. Such structures are known as Galton boards 
or Lorentz gas. According to the mathematical results 
of Sinai the dynamics on such a lattice is completely 
chaotic 0|- The theoretical studies ^( performed to 
understand these experiments showed that the chaotic 
classical dynamics and periodic orbits significantly affect 
the conduction properties in such superlatticies. 

The effects of microwave radiation on the conduction 
properties of antidot superlattices has been addressed in 
experiments jlfij . However in these structures due to the 
symmetry of the superlattice the ratchet effect was for- 
bidden. Asymmetric mesoscopic structures under exter- 
nal periodic driving have been addressed in the experi- 
ments 0| • However zero mean force ratchet was absent 
due to low frequency of driving which was essentially adi- 



abatic (lif . 

The recent theoretical studies of dissipative transport 
in asymmetric superlatticies showed that microwave ra- 
diation induces directed transport in such systems (zero 
mean force ratchet) [l9lEoT|. These works were mainly 
performed for a semidisk Galton board which is obtained 
from the usual Galton board of rigid disks by replacing 
each disk with semidisk oriented in one fixed direction 
(see Fig. 1). In Ref. the model of a friction force 
tf = — 7717V with constant friction coefficient 7 has been 
used for a particle of mass m moving with velocity v. In 
Ref. [2(J the case of particles in a Maxwell thermostat 
at temperature T was considered. It was shown that 
the thermostat creates an effective friction coefficient 7 
which depends on the microwave field strength and the 
temperature of the thermostat. However the most in- 
teresting case is the Fermi-Dirac thermostat since it de- 
scribes the experimental conditions of antidot superlat- 
ticies with 2DES |2l|. Until now no numerical studies 
were performed in this regime. Only theoretical esti- 
mates have been proposed in Ref. |20j. Their validity 
was never checked and remains questionable. 

In this work I develop a numerical method which allows 
to study the directed transport induced by polarized mi- 
crowave radiation f = /(cos 9, sm9) cos(cut) in 2DES at 
various values of the Fermi energy Ep and temperature 
T (here w is the radiation frequency and 9 is the polar- 
ization angle with respect to the x axis in Fig. 1). On 
the basis of this method I performed extensive numeri- 
cal studies which allowed to establish the dependence of 
ratchet flow velocity Vf on various system parameters in- 
cluding T and Ep. Contrary to the estimates proposed 
in Ref. [2£j the dependence on T is weak when T <C Ep. 
The obtained results allow to predict typical values of 
currents in asymmetric antidot superlattices. 

The paper is organized as follows, in Section II the 
description of the model and of the numerical method is 
presented, the results are described in Section III, and 
conclusion is given in the last Section. 
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FIG. 1: (color online) Example of the semidisk Galton board 
with one chaotic trajectory in the region —5.5 < x < 1.5 and 
19 < y < 26. Semidisk scatterers, separated by a distance 
R = 2, are shown in black. The trajectory is displayed in red 
/ gray. The strength of microwave radiation and its frequency 
are / = 5, u = 1, polarization angle 8 = and temperature 
T = 0.1-Ef. Here and in other figures Ef = 10 and disk radius 
ra = 1- The Metropolis algorithm parameters are AE/Ef = 
0.075, At = 0.005. 
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FIG. 2: (color online) Energy distribution functions fn(E) ob- 
tained numerically for different field amplitude values / and 
different temperatures T. For the green, black and red curves 
the temperatures are T/Ef = 0.01, 0.1, and 0.4 respectively 
(at E I 'Ef = 1.1 the order of curves is from bottom to up). 
The applied force / is zero for the dashed curves and / = 5 for 
the full curves (u = 1, = 0, R = 2). At / = the dashed 
curves coincide with the Fermi-Dirac distribution /f(-B) at 
corresponding temperature T. The Metropolis algorithm pa- 
rameters are the same as in Fig. 1. 



MODEL DESCRIPTION 

The geometry of the model is represented by a Galton 
board of semidisks which form a two-dimensional hexag- 
onal lattice as represented in Fig. 1. The radius of the 
semidisk is r d and the distance between the disk centers 
is R. A particle with mass m moves under the action of a 
force f created by a polarized microwave radiation. The 
collisions with the semidisks are elastic. In the numerical 
simulations I use dimensionlcss units with m = r d = 1 
and the Fermi energy Ep — 10 (thus the Fermi velocity 
V F = v20) • In order to convert the numerical results ob- 
tained with this units, one should determine the dimen- 
sionless ratios: for example in Fig. 1 the frequency u) = 1 
and field strength / = 5 corresponds to dimensionless 
values uj = ur d /V F = I/a/20- and / = fr d /E F = 1/4. 

It is assumed that particles are non interacting but 
that the contact with a thermostat creates the Fermi- 
Dirac distribution f F (E) = l/(E F [exp((E-E F )/T) + l}) 
at temperature T and Fermi energy Ep, where E is 
a particle kinetic energy. Here the energy distribution 
for one particle f F {E) is normalized by the condition 
J fp(E)dE — 1. This one particle distribution gives 
also the result for many particles by simple rescaling by 



the number of particles as it is usually done for 2DES 
(see [2l| p. 193). Such a situation corresponds to exper- 
iments with 2DES in antidot lattices similar to those of 
Refs. To keep particles in a thermal equilib- 

rium with the Maxwell distribution it is possible to use 
various methods including the Nose-Hoover thermostat 
used in Ref. [2(|. However for the Fermi-Dirac thermal 
distribution this method is not appropriate and a new 
approach should be developed. Indeed the Nose-Hoover 
equations are constructed in such a way that th ey g ive the 
Maxwell distribution in particle velocities [2(], l22| . They 
should be significantly modified to generate the Fermi- 
Dirac distribution and until now this problem has not 
been addressed yet. The first attempt by the authors of 
Ref. [2(J did not succeed in achieving a good convergence 
to the Fermi-Dirac distribution |2^|. 

My approach is inspired by the successful Monte Carlo 
method adapted to simulate numerically the transport 
properties in semiconductor devices 24] . To obtain a 
stable Fermi-Dirac thermalization of particles on the 
semidisk Galton board (see Fig. ^) the following pro- 
cedure has been applied: (i) the equations of motion 
were integrated exactly on the time interval At using 
the analytical solution; (ii) after that the particle en- 
ergy is changed from its value E to another value in the 
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FIG. 3: (color online) Directed transport for one trajectory at 
various polarizations of radiation 8 = 0, tt/5, 7r/4, n/2 (from 
left to right clockwise), other parameters are set as in Fig. 1. 

interval (E — AE, E + AE) without changing the direc- 
tion of the particle momentum. The choice of this value 
is governed by the Metropolis algorithm [2^] which im- 
poses the convergence to the Fermi-Dirac distribution fp. 
Namely, a random value E' is chosen in (E — AE, E + 
AE), with probability mm(f F (E')/f F (E),l) the algo- 
rithm sets E = E' , otherwise E remains unchanged. Af- 
terward the algorithm returns to step (i). The times 
when collisions with semidisks occur are found with the 
precise Newton algorithm as in Ref. 0. The step AE 
can be considered as a thermalization step which deter- 
mines the rate of convergence to the equilibrium distri- 
bution fp. 

The Metropolis algorithm described above gives con- 
vergence to the equilibrium distribution fp in the ab- 
sence of microwave radiation. The examples of steady 
state distributions at different temperatures T are shown 
in Fig. 2. The proximity of the numerically obtained 
distribution f n (E) to the theoretical steady state /f(E) 
can be characterized by the dimensionless mean square 
deviation a = E F J(f„(E) - f F {E)) 2 dE. This quantity 
remained small in all numerical simulations at / = 
showing a good convergence to the Fermi-Dirac distribu- 
tion (e.g for the cases of Fig. 2: a = 1.8 x 1(T 5 , 2.3 x 1(T 5 
and 8.3 x 1(T 5 for T/E F = 0.4,0.1,0.01 respectively). 

However, the introduction of the microwave radiation 
modifies the distribution f n which depends on / and 
other system parameters. This is clearly seen in the typi- 
cal cases presented in Fig. 2. For relatively high temper- 
atures the distribution /„ remains a smooth monotonic 
function of energy whereas at low T the microwave field 
creates a characteristic peak near the Fermi energy Ep. 
As a result the developed numerical method allows to 
study the transport created by microwave radiation in 
an asymmetric semidisk Galton board in the stationary 
regime. To be close to realistic experimental situations 
additional random scattering has been introduced to take 
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FIG. 4: (color online) The ratio of flow velocity Vf to the 
Fermi velocity Vf as a function of the thermalization step 
AE for a fixed At = 0.005. The different curves correspond 
to values of /: / = 7.0, 6.0, 5.0, 4.0, 3.0 with respective colors: 
violet, blue, green, red, black (from top to bottom at AE/Ef 
= 0.2). The symbols mark the temperature T/Ef = 0.01 
(squares) and T/Ef =0.1 (circles). Other parameters are as 
in Fig. 1. 



into account the effect of impurities. Namely after time Tj 
the direction of particle momentum is changed randomly 
(angle changes in the interval [0, 2tt]). In the majority 
of cases studied the value Tj was kept sufficiently large 
{jiV F /rd > 1000) and did not influence the stationary 
transport properties (the dependence on n will be dis- 
cussed in the next Section). Here V F = \ / 2Ep is the 
Fermi velocity. 

This stationary regime, which sets in presence of ra- 
diation, clearly shows the photogalvanic (ratchet) effect 
with directed transport born from chaos as it is shown 
in Fig. 3 for different polarization angles 9. I would like 
to stress that this effect originates from the bulk of the 
sample; it is known that radiation can lead to the appear- 
ance of electron-hole pairs, which may diffuse to different 
contacts because of the electric field in the depletion re- 
gion of the junction, however in this case the effect is 
due to the contacts and not to the bulk of the sample. 
For 6 = the transport is directed to the left (negative 
direction on the x axis), while for = tt/2 the transport 
is oriented to the right. The polarization dependence is 
similar to that obtained in previous works |13 . |2fJ and is 
well described by the relation ip w it — 29 where ip is the 
angle between the transport direction and the x axis. In 
this way the average velocity of transport can be written 
as Vf = Vf(cosip, sin ip). It is calculated by following a 
single trajectory for a long typical time Wp/rd ~ 10 7 . 
It was also checked that the averaging over an ensemble 
of few tens of different trajectories gives statistically the 
same result. 

Special checks have been made to ensure that the 
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ratchet velocity v / does not depend on the Metropolis al- 
gorithm parameters. An example of such checks is shown 
in Fig. 4 where the thermalization step AE was changed 
by one order of magnitude. In spite of this variation 
the value of Vf remains stable. It was also checked that 
the resulting ratchet velocity is not sensitive to varia- 
tion of At, e.g a variation of At by more than one order 
of magnitude gave no variation of Vf within a 5% accu- 
racy. Physically it is possible to say that the relaxation 
time to the equilibrium r re i is approximately given by 
the relation r re i ~ AtEpj AE 2 . Thus, the numerical 
checks above can be physically interpreted as the inde- 
pendence of the ratchet velocity on the variation of the 
energy relaxation timescale T re i which varied by two or- 
ders of magnitude. A similar effect has been seen with 
the Nose- Hoover thermostat for the Maxwell equilibrium 
distribution [2(| . This result is also in a qualitative argu- 
ment with the theoretical arguments given in Ref. 0, L3] 
according to which r re i does not directly affect Vf. In my 
further simulations the Metropolis algorithm parameters 
are set to typical values AE/E F = 0.075, At = 0.005. 
The data shown in Fig. 4 clearly demonstrate that v/ 
grows with the radiation strength /. At the same time 
Vf is not very sensitive to the variation of temperature 
T. Detailed studies of parameter dependence of Vf are 
presented in the next Section. 

The results presented above show that the developed 
algorithm allows to simulate non interacting electrons 
at thermal equilibrium with the Fermi-Dirac distribu- 
tion. The microwave driving is assumed to be relatively 
weak so that it gives only small deviations from the un- 
perturbed distribution, this fact is clearly illustrated in 
Fig. 2. In this regime the perturbed distribution of parti- 
cles is determined by the unperturbed distribution fp (E) 
and the microwave field, the effect of the latter is exactly 
taken into account by the Hamiltonian equations of mo- 
tion. Therefore the developed algorithm should correctly 
describe the non equilibrium steady state distribution 
that emerges under mircrowave driving. This is indi- 
rectly confirmed by the fact that the directed transport 
is not sensitive to the thermalization step of the Metropo- 
lis algorithm (see Fig. 4). In a sense the Metropolis steps 
combined with the Hamiltonian equations of motion give 
the solution of the kinetic Boltzmann equation in the 
presence of microwave driving. 

The above arguments should be also valid for another 
unperturbed thermal distribution, e.g the Maxwell dis- 
tribution fm(E) = exp(—E/T)/T. This case was an- 
alyzed in Ref. j20J on the basis of Nose-Hoover equa- 
tions. In fact the Metropolis algorithm had been invented 
to treat the Maxwell thermal equilibrium J2{|. Thus, I 
made numerical tests with the Metropolis algorithm for 
the Maxwell distribution. The obtained results repro- 
duce the functional dependences found in Ref. [13] (see 
Eqs. (4,5) there) with approximately the same values of 
the numerical constants. This gives independent confir- 



mation that the Metropolis algorithm treats correctly a 
weak external perturbation that drives the system out of 
equilibrium. It also shows that various thermal distribu- 
tions can be treated by this method. 

The model I described assumes that the electrons are 
non interacting. To be valid it requires that the Coulomb 
interaction between electrons in 2DES E ee « e 2 ^/nn e /e r 
is small compared to the kinetic energy given by Ep = 
irn e H 2 /m. Here n e is the electron density, e r is the di- 
electric constant, e is the electron charge and m is the 
effective electron mass. Thus the effective strength of 
interaction is characterized by the dimcnsionlcss param- 
eter r s = E ee /Ep |26( which should be small. Its value 
for experimental 2DES obtained in GaAs/AlGaAs het- 
erostructure with electron densities n e sa 10 12 cm -2 , an 
effective electron mass m ~ 0.065m e and dielectric con- 
stant e r = 13 is approximately r s ~ 1. At such values 
the interaction between quasiparticles is considered to 
be weak and is usually neglected in a first approximation 
[2lH26| . for example the Wigner crystal typically appears 
at r s i=s 37. 

At such densities n e inside a cell of size S = Ifim 2 
the quantum level number of an electron at the Fermi 
energy is Np = n e S ~ 10 4 . Therefore electrons are in a 
deep semiclassical regime and the classical Monte Carlo 
approach used above is well justified, see also |24| . 

NUMERICAL RESULTS 

I have investigated the dependence of Vf on several 
system parameters which are relevant for a realistic ex- 
periments with 2DES in antidot lattices. Among them 
are the temperature T, the microwave field strength / 
and the microwave frequency ui. The effects of geome- 
try are studied by changing the lattice constant R that 
allows to choose the optimal regime where the photogal- 
vanic effect is stronger. The effects of impurity scattering 
is modeled by variation of scattering time that gives 
insight on the stability of the effect in respect to experi- 
mental imperfections. At last the effect of magnetic field 
is also analyzed. 

The temperature dependence for a typical set of pa- 
rameters is given in Fig. 5. The obtained numerical data 
show that there is a weak drop of the ratchet velocity 
Vf with the increase of temperature T. However in the 
regime with T <C Ep, Vf is practically temperature in- 
dependent. This dependence is preserved for various ra- 
diation strengths /. The velocity of transport increases 
with the growth of /. A detailed study of the effect de- 
pendence on the microwave field / is presented in Fig. 
6. It shows that the dependence on / is quadratic in the 
region T <C Ep, and temperature independent in a large 
interval of field strength. At higher / a deviation from 
the quadratic dependence starts to be visible, this devi- 
ation starts earlier at high temperatures. Thus on the 
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FIG. 5: (color online) Top panel: dependence of the rescaled 
flow velocity Vf/VF on the rescaled temperature T/Ef, dif- 
ferent curves correspond to / = 7.0, 6.0, 5.0, 4.0, 3.0 (from 
top to bottom). Bottom panel: the same data are shown as 
a function of E F /T. Here u) = 1, R = 2, 6 = 0. 

basis of the results presented in Figs. 5,6 it is possible 
to conclude that at T <§; E F and in the limit of weak 
driving the flow velocity is described by the relation : 

v f /V F =C(r d f/E F ) 2 (1) 

Here the dimensionless factor C may depend on the mi- 
crowave frequency, lattice geometry, and impurity scat- 
tering time. However it is independent of T and /. For 
R = 2, LorijVp -C 1 and nVp/rd 3> 1 the obtained data 
give C = 0.129 ±0.002. 

The dependence (1) is qualitatively different from the 
theoretical estimates proposed in Ref. [2(|- To under- 
stand the origin of this difference I remind the main el- 
ements of estimates given in Refs. 0, H^- They are 
based on the fact that the microwave radiation produces 
diffusive energy growth of electron energy in time with 
the rate : Dp = (5E) 2 /5t. Here SE is the energy vari- 
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FIG. 6: The rescaled flow velocity Vf/VF as a function 
of rescaled applied force rdf/EF for different temperatures: 
T/Ef = 0.4 (green diamonds), 0.1 (blue squares), 0.01 (black 
circles). The red (full gray) curve shows a parabolic fit of 
data vf/VF = C (r d f/E F f with C = 0.129 ± 0.002 at 
T/Ef = 0.01 (the fit is done in the interval [0,0.45]). Here 
u = 1, R = 2, e = 0. 



ation after a time St. In the limit of low frequency 
driving it is possible to write Dp ~ (E) 2 t c ~ f 2 V F l 
where E is the particle energy and I is the mean free 
path which is I ~ R 2 /r<j. If a particle would experi- 
ence a friction force f/ = —7717V, the diffusion in ene rgy 
would give energy variation (SE) 2 ~ Dp/j. In Ref. |20j 
it was assumed that SE is fixed by the thermal distri- 
bution so that SE ~ T. Thus the statistical station- 
ary distribution imposes an effective friction with coef- 
ficient 7 ~ Dp/T 2 . This relation is important because 
the ratchet velocity is given by relation Vf ~ Z7 in the 
limit of weak friction as it has been shown in Ref. 0] by 
extensive numerical simulations. For the Maxwell distri- 
bution this gives Vf ~ (lf) 2 /T 3 ^ 2 m 1 ^ 2 since in this case 
Dp ~ / 2 (T/m) 1 / 2 Z, where it is used that the thermal ve- 
locity v ~ s/T/m. The extension of the above arguments 
to the Fermi-Dirac distribution led the authors of to 
the conclusion that vj/V F ~ C(fl/T) 2 . However this 
expression is drastically different from Eq. (1) according 
to which there is no temperature dependence in the limit 
T <C E F . This contradiction can be resolved if the varia- 
tion of particle energy induced by the radiation is of the 
order SE ~ E F (and not SE ~ T). Indeed in the free 
electron model the particle can move in the whole energy 
interval defined by the Fermi energy whereas in the esti- 
mates in Ref. |20| it was assumed that the particle can 
move only in the narrow thermal layer near the Fermi 
surface. In the free particle model it is therefore rather 
natural that the particle energy variation is SE ~ E F , 
that leads to the result of Eq. (1). 
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The fact that the free electron model remains valid 
in the presence of microwave driving can be also under- 
stood from the following arguments. For non interacting 
electrons the Hamiltonian is given by the sum of one 
particle operators (microwave driving is also one particle 
operator). Hence, the many particle state (wave function 
or density matrix) is obtained simply from one particle 
states by antisymmetrization. Thus, the Pauli princi- 
ple can be taken into account by averaging the final one 
particle results over the Fermi-Dirac distribution. This 
statement also explains the validity of the classical ki- 
netic Boltzmann equation for the description of transport 
properties of metals. It is demonstrated more rigorously 
for a non interacting Fermi gas in (Chap. 5 in Ref. |2l| 
and Ref. [13|)- As a consequence the Pauli blockade does 
not appear for non interacting particles and the argu- 
ments presented in Ref. [2(| are not valid at least for 
weak r s values. 

The arguments given above allow to understand the 
physical origin of the relation between the ratchet veloc- 
ity and effective friction coefficient induced by microwave 
radiation that perturbs the system out of thermal equi- 
librium. Another approach has been developed recently 
in Ref. |28J using perturbation theory for the Boltzmann 
kinetic equation in the limit of weak radiation and weak 
density of randomly distributed asymmetric scatterers. 
The model proposed in Ref. [2^] is rather different from 
the one considered here, e.g. scatterers are distributed 
randomly on the plane, their density is required to be 
small and impurity scattering is necessary for the regu- 
larization of the model. However in spite of these differ- 
ences the global dependence of ratchet velocity on Fermi 
energy Ep and radiation strength / is the same as in Eq. 

(I)- ' 

The frequency dependence of vt /Vf is shown in Fig. 7. 
The data presented there demonstrate that the frequency 
spectrum is independent of temperature (for T <C Ep) 
and that quadratic dependence on / is valid for a large 
frequency range. The spectral dependence has a few 
characteristic features. For ujrd/Vp -C 1 the ratchet ve- 
locity becomes independent of u, that is in agreement 
with the fact that Dp oc Vf is independent of u>. This 
has been also seen in models analyzed in Refs. [TflL Enj|. 
For uird/Vp ^> 1 the ratchet velocity drops with u>. This 
is consistent with the dependence of Dp on w which in 
this limit can be estimated as Dp ~ f 2 Vp/uj 2 l. This 
comes from the fact that at high frequency the change of 
velocity after one collision with semidisks is Av ~ / /com 
so that the energy is changed by 5E ~ Vpf /u> after a 
time between two collisions 5t ~ l/Vp (I ~ for R ~ 
and I ~ R 2 /rd for R ^> r^). The frequency dependence 
shown in Fig. 7 has also a resonance at uird/Vp » 3. This 
is probably related to high frequency collisions in narrow 
bottle necks between two semidisks (see Fig. 1 where the 
narrow bottle neck is approximately by a factor of three 




FIG. 7: (color online) Dependence of rescaled flow veloc- 
ity Vf/VF on the rescaled microwave frequency wrd/Vp, for 
T/E F = 0.01, / = 5.0 (red squares); T/E F = 0.1, / = 5.0 
(green circles); and T/Ef = 0.01, / = 3.0 (black triangles, 
in this case Vf is multiplied by factor (5/3) 2 to underline 
quadratic dependence on /). Here R = 2, — 0. 
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FIG. 8: Dependence of V//Vf on the rescaled distance be- 
tween semidisk centers R/rd for different frequency values 
U) = 1, 1.5, 2.0 (curves from top to bottom). The smooth 
curve shows the fit of Eq.(2) for u) = 1 with A = 0.017, 
B = 0.012. Here / = 5.0, w = 1, 9 = 0, T/E F = 0.1. 



smaller than disk radius). 

The dependence of Vf on the distance between disk 
centers R is shown in Fig. 8. Initially Vf starts to grow 
with R then reaches a maximum value and drops at large 
R values. The position of the maximum depends on the 
microwave frequency. With the increase of frequency the 
maximum moves to smaller values of R. Qualitatively 
this corresponds to the situation when the microwave fre- 
quency becomes comparable with the frequency of colli- 
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sions of particles with semidisks. The dependence of data 
on parameters can be satisfactory described by a fit for- 
mula : 



v f /V F = A 



2 f2 



R f 



1 



E 2 

F 



1 + B{LoR 2 /r d V F ) 2 



(2) 



Here A, B are dimensionless fitting parameters. The fit 
for three values of u in Fig. 8 gives A i=s 0.017 and 
B ps 0.012. For u> — ► this expression is in a satisfac- 
tory agreement with the value C « 0.13 found in Fig. 6 
at R = 2. The physical origin of this fit is related to 
frequency dependence of the diffusion rate De which in- 
terpolates between the low frequency regime {De inde- 
pendent of to) and the high frequency regime where De 
drops quadratically with u> (see estimates given above). 
Eq. (2) gives reasonable description of obtained numeri- 
cal data in the regime where I is not too large compared 
to R. This situation is most interesting for direct exper- 
imental studies where R is not very large compared to 

Another important experimental parameter is the scat- 
tering time induced by impurities which are always 
present in real samples. The data presented in previ- 
ous figures were obtained in the regime of very large 
r,-. The effect of finite values of Tj on the ratchet ve- 
locity is described in Fig. 9 for various lattice constants 
R. The data show that at large n values Vf is indepen- 
dent of the impurity scattering time while at small Tj, 
vj/V F drops approximately linearly with T^Vpfrd- In- 
deed the asymmetry of semidisks is washed out by im- 
purity scattering and the ratchet effect should disappear 
at small 77. At the same time it is important to stress 
that the presence of impurities is not necessary for the 
onset of directed transport. In experimental conditions Tj 
can depend on temperature because of electron-phonon 
scattering or electron-electron interactions which give a 
dependence of ti on T . This may lead to a significant 
temperature dependence of the photogalvanic effect, in 
the temperature range T ~ 1°K . 

Experimentally it is also possible to study the depen- 
dence of the effect on magnetic field B perpendicular 
to the 2DES plane. To investigate this dependence the 
method described above was adapted to the presence of 
a magnetic field, which was included in the analytical 
solution of motion equations between Metropolis ther- 
malization steps. The magnetic field dependence is given 
in Fig. 10. The data clearly shows that the ratchet effect 
disappears for sufficiently strong magnetic fields when 
the Larmor radius r; = V F / B becomes smaller than the 
distance R between semidisks (here electron charge and 
mass are set to 1). Indeed for 77 <C R the classical elec- 
tron dynamics becomes integrable and the diffusion rate 
in energy De goes to zero due to absence of chaos, thus 
leading to the disappearance of ratchet {vf oc De)- In 
principle the magnetic field changes the transport direc- 
tion (angle ip). I do not discuss this dependence here 




FIG. 9: (color online) Dependence of Vf/VF on the rescaled 
impurity scattering time VFn/rd for R/ra = 6.0, 4.0, 3.0, 
2.0 (from top to bottom respectively at VFTi/ra = 50). Here 
/ = 5.0, w = 1, e = 0, T/E F = 0.1. 



0.08 



0.06 



> 



. 0.04 



0.02 




r d B/V F 

FIG. 10: (color online) Dependence of Vf /Vf on the rescaled 
magnetic field r^B/VF for R/rd = 4.0, 3.0, 2.0 (curves from 
top to bottom at r d B/V F = 0). Here / = 5.0, u = 1, = 0, 
T/Ef = 0.1. 



since the main point is that the ratchet effect disappears 
at relatively low B (see below). 



CONCLUSIONS 

The obtained results clearly shows the existence of zero 
mean force ratchet in asymmetric semiconductor struc- 
tures induced by microwave radiation. The obtained re- 
sults give the following dependence for the strength of 
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the stationary current induced by the ratchet effect in 
one row of semidisks (row width \^3R) : 

T R R A ^ emR3 m 

I = V3en e Rv f = A^j — -^ (3) 

where Ep = Trn e h 2 /m and m — 0.065m e . This depen- 
dence holds in the low frequency regime which is usu- 
ally satisfied at typical electron densities n e = 10 12 cm -2 
where the collision frequency is of the order of 200GHz 
for a R ~ 1/zm and under the assumption R ~ r<i For 
R = Ifim and //e = lV/cm, the equation (3) gives the 
current O.lnA. In samples with high mobility the mean 
free path can have values as high as 5/im, and therefore 
the optimal regime for photogalvanic effect will be when 
R is of the same order and it is quite possible that in this 
situation the current per row can be as high as lOnA. Ac- 
cording to the results of Fig. 10 for R ~ 1/mi the ratchet 
effect starts to disappear at magnetic field B ~ 0.1T. 

The asymmetric antidot lattice can be considered as a 
prototype for transport in asymmetric molecular struc- 
tures. The latter have attracted recently a significant in- 
terest in view of possible biological applications of ratch- 
ets [i^. Therefore experimental investigations on the 
ratchet effect discussed in this paper are highly desirable. 

I thank D. L. Shepelyansky for stimulating discussions 
and for his interest in this research. I am also grateful 
to M. V. Entin, and L. I. Magarill for access to their 
unpublished results on their kinetic equation approach 
to the photogalvanic effect. This research is done in the 
frame of the ANR PNANO project MICONANO. 
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